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ABSTRACT 


This research indicates that a knowledge of angle of attack (o) 
and/or gust perturbation angle of attack ία»; is of prime importance 
for a longitudinal autopilot designed to alleviate the effect of 
vertical gusts. However, accurate measurement of о апа = is extremely 
difficult. This research indicates that the readily measurable normal 
acceleration factor (n) and pitch rate (q) can be used to provide a 


and = information indirectly. 
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SYMBOLS 
Linear acceleration of aircraft center of mass in 
z direction, ft/sec. 


Length of mean aerodynamic chord, ft. 


0C 
m 


on С 


iy 

δα ς 

local acceleration of gravity, ft/sec. 
μυς, {εξ 

Mean-square value of output of system. 
Control system gains. 

Gust characteristic length, ft. 
Aircraft mass, Slugs. 


a. /g, normal acceleration factor, ft/sec?/g. 
2 


He eae = n2 , (ft/sec?/g)?. 


Value of n? with no autopilot utilized, (ft/sec?/g)?. 
Pitch rate, rad/airsec. 

Pitch rate due to vertical gust, rad/airsec. 
Aircraft forward velocity, ft/sec. 

Aircraft forward velocity perturbation, ft/sec. 
Stability axis system. 

Velocity perturbation along z axis, ft/sec. 

Gust vertical velocity, ft/sec. 


WA angle of attack perturbation due to gust, 


measured at aircraft center ОЁ mass, rad. 





Elevator deflection, measured from trim, rad. 


- Mean-square value of elevator deflection (Rad)?. 


DIE (direct lift control) deflection. 
Throttle deflection. 
Spatial frequency, rad/ft. 


QU t* = QU. = , Frequency, Rad/AirSec. 
O 


Root mean-square value of vertical gust velocity, ft/sec. 
Pitch angle perturbation, Rad. 


Atmospheric density, Slugs/Ft?. 
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I. INTRODUCTION 


A. LONGITUDINAL AUTOPILOTS 
1. Basic Ideas in Longitudinal Control 

Since the beginning of powered flight, methods have been sought 
to reduce pilot effort in controlling aircraft and to automate various 
aspects of aircraft control. The autopilots and stability augmentation 
systems of today are realizations of such research. 

Autopilots and stability augmentation systems are usually 
classified as longitudinal or lateral. As the names imply, longitudinal 
(lateral) autopilots are concerned with the longitudinal (lateral) flight 
regulation. This dual classification is justified in flight conditions 
where the existence of pure longitudinal/lateral motion is ensured. 

Blakelock [Ref. 1] offers a very readable introduction to auto- 
pilot design. He discusses three basic types of longitudinal controllers, 
each designed for different longitudinal operations: a displacement 
autopilot where aircraft pitch angle is fed back and compared with a 
reference pitch angle; a pitch orientational autopilot where aircraft 
pitch rate is fed back and compared with some desired pitch rate; and an 
acceleration control system where the acceleration of the center of 
gravity is fed back and compared with a reference acceleration. A good 
treatment of the feedbacks that may be used in longitudinal control 
systems can be found in Ref. 2, Sec. II. Table I (obtained from Ref. 2) 
indicates the different feedbacks which can be utilized to control 


longitudinal motion. 





The research summarized as is concerned with autopilots designed 
to alleviate the effect of vertical gusts on a large jet transport in 
steady, wings-level flight. Two flight conditions are considered, high 
altitude cruise and landing approach. By reducing the aircraft's gust 
response it is possible to provide for a smoother ride, aircraft 
structural fatigue problems can be alleviated, and the ability of the 
aircraft to function as a gun platform can be imporved. 

2. Iongitudinal Control to Alleviate Gusts 

Appendix B offers an introduction to the statistical design 
concepts utilized in this study. As Table I indicates many feedback 
parameters such as pitch, pitch rate, center of gravity acceleration, 
etc., can be utilized as feedbacks in a controller designed to reduce 
gust response. Hess [Ref. 3] discussed a two-parameter gust alleviation 
system shown in Fig. l. He analyzed this system using a performance 


index given by: 


= r=i So q | 一 (s) | 2de + Se, D Ka E 
co "o S=Ju E S=Jw 


where В is the non-dimensionalized frequency and the transfer functions 

a; 9) and SEN are obtained fram the aircraft's equation of motion as 
discussed m Appendix C. d 15 the non-dimensionalized power spectral 
density of the gust field and was obtained from Fguation 1 found in 
Appendix B. Hess's study utilized a large, rigid jet transport. The 
characteristics for the two flight conditions considered may be found in 


Appendix A. His results are shown in Table II. 


10 





Hess also proposed an optimal control system, Fig. 2, which 
showed considerable improvement over the two-parameter feedback system 
above. For analysis of this system Hess used a performance index given 
by: 


_ Lim 
E [qn? (t) + rut (t)]at 


HI 
oh, 


The aircraft was the same as used in analysis of the two-parameter system. 
Table III shows results of his study for a a/r equal to 10, a value was 
selected by a root-square locus technique. Note the improvement in 


An? which Hess's optimal system maintains over the two- 


D^unalleviated 


parameter system of Table II. 


B. BASIS OF THIS RESEARCH 

A major obstacle to the implementation of Hess's optimal gust 
alleviation system is the difficulty in measuring Sec Attempts to design 
instruments that will give satisfactory on-line measurements Of SÉ have 
been unsuccessful. This report investigated the effect of eliminating 
= information on the performance of the alleviation system. Figure 3 
shows Hess's optimal gust alleviation system with gust information 


removed. This is the system utilized in this study. 


11 





II. SOLUTION PROCEDURES 


A. METHODS OF SOLUTION 

Optimal feedback design can be classified as free configuration or 
fixed configuration. The latter is often referred to as parameter opti- 
mization. In the free configuration procedure the designer is not 
restricted in any way. The final design depends only upon the input and 
the form of the performance criteria. With the parameter optimization 
techniques the configuration is chosen using feedbacks which are readily 
obtainable and then the feedback parameters are optimized to minimize 
the performance criteria. 

As noted in Ref. 2, for autopilot design the parameter optimization 
technique is probably the superior method considering current state of 
the art. This technique yields an autopilot which fits the "real world" 
situation. This research utilized the parameter optimization procedure 


since the form of the longitudinal controller was prespecified (Fig.-3). 


B. PERFORMANCE INDEX 
The performance index chosen for the optimization procedure was 


given by: 
WA Cn C A USA A 
I= E + 5 EN = n“+ n (1) 
where ш denotes the non-dimensionalized frequency of the randam gusts. 


This performance index is identical to that used by Hess in his evaluation 


of the optimal controller with q/r = 10. It is a sum of the mean square 
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normal acceleration factor and the mean square elevator deflection. The 
latter was included to prevent large elevator excursions. Equation (1) 
can be manipulated as indicated in Appendix B to a more usable form: 


со 


r=} Se, QU) Eom L Se, 08) 12- "Deia (2) 


со 


Sa is obtained from Equation (1) in Appendix B by substituting w = qu t*. 


"his non-dimensionalized form is: 


2 


a 








Ф e κα c : 1 (3) 
og Léo. e? (©?/412) 


The transfer functions 2 (ш) апа ` (ju) are obtained from the air- 
craft's longitudinal rs ee the equation for normal 
acceleration and the mathematical equation of the autopilot. Appendix C 
shows how this was done. 

For the purposes of determining the aerodynamic gust effects, the 
aircraft was idealized as a segment of the x-axis. This idealization 
has been found to give good results as long as the shortest vertical- 
gust wave length is greater or equal to eight times the aircraft length 


Mer. 4, pg. 323]. 


C. INTEGRATION OF THE PERFORMANCE INDEX 

Substitution of Equation (3) and the transfer functions into 
Equation (2) produced an equation whose form is frequently found in system 
analysis. Tables for solution of these equations are readily avallable, 


era, Ref, >. 
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Appendix D shows how a computer program was written to solve for 
the performance index, which in this case took on the functional form 


Of Ig (K1,K2,K3,1, )- 


р. OPTIMIZATION 

Integration of the performance index for a given characteristic gust 
length (L produced a function which contained the variables Kl, K2, 
and K3, i.e. the feedback gains. The function obtained could be thought 
of as a performance index surface with dimensions Kl, K2 and K3. The 
optimm solution would be found where the surface produced a minimum 
value. A computer program was written to solve for this minimum by a 


method of steepest descent. Appendix D indicates how this was accomplished. 
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III. RESULTS AND DISCUSSTON 


A. PERFORMANCE COMPARISONS 

The results of the computer computation are presented in Table IV. 
The performance surface at most points was very shallow making it diffi- 
cult to rapidly find a minimum value on the surface. It was found that 
the performance index was a weak function of K2, but a strong function 
ӨТ КЇ, K3 and Ly Kl, K2 and K3 were found to be very dependent on the 
flight condition, i.e. cruise or landing approach. 

Comparing Table IV with Tables II and III shows that the three- 
parameter system has the largest perfomance index. Figures 5 and 6 


indicate that the An*/n2 performance initially compares 


unallev. 
favorably with the optimal and the two-parameter system, but the per- 
formance rapidly deteriorates as L, increases. This indicates that a 
knowledge of = becomes increasingly important as L. gets larger. 

Note that in the two-parameter system of Fig. l the feedbacks n and 
q give a knowledge of a and SC This Can be seen Iram the eguation for 


the normal acceleration factor which can be written as: 


D(a + a ) = q - 
(а Zei q 


NIA 
Q 
A 
名 


B. TRANSIENT PERFORMANCE 

Figures 6 and 7 present the locus of roots of the characteristic 
equation obtained when the feedback gains Kl and K3 are set at values 
obtained for a gust characteristic length of 1000 feet while K2 was 


allowed to vary. This could be done because the performance index was 
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a weak function of K2. A value of the variable K2 could be selected 
which would give the three-parameter system better transient performance 
than the two-parameter system, but this transient performance was still 
inferior to the optimal control system. A new root locus diagram would 
have to be drawn for each gust length and flight condition. In an effort 
to circumvent the drawing of multiple root locus diagrams, the feedback 
gains Kl and K3 were set at values obtained for L. equal to 3000 feet. 

It is found, however, that by doing this performance was unacceptably 


degraded, particularly in the lower ranges of Lo 


C. STABILITY 
It was found that the three-parameter system was stable throughout 


the range of feedback gains considered here. 


D. POSSIBLE ALTERNATE SOLUTION 

Comparing the three-parameter longitudinal controller of this 
research with the optimal controller (Fig. 2) and the two-parameter 
System (Fig. 1) it appears that I is an important parameter for a gust 
alleviation system. Since, using present technology, it is not possible 
to measure = aboard the aircraft some alternative method of obtaining 
this information must be employed. 

It was noted previously (Sec. III, Part A) that a knowledge of n and 
a would provide information on a and x. This would explain why the RS 
performance of the two parameter feed back system tends to approach the 
RMS performance of the optimal controller as the characteristic gust 


length increases. 
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thy 


It is felt that by combining the optimal controller and the two- 
parameter system a gust alleviation system could be built that would have 
performance comparable to the optimal controller. Indeed, Hess [Ref. 6] 
has shown this to be the case. The modified gust alleviation system 
appears in Fig. 4. This system would have the advantage of utilizing 


feedback information that is readily obtainable on-board the aircraft. 
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IV. CONCLUSIONS 


l. The longitudinal gust alleviation system evaluated by this report is 
a stable, usable system and could be designed to provide better transient 
performance than the two-parameter system of Fig. 1. However, the two- 
parameter system and the optimal controller of Fig. 2 both exhibit better 
RMS performance. 

Dr E provides important information for a gust alleviation system. 
Knowledge of > becomes increasingly important as the gust characteristic 
length increases, i.e. as the random gust field departs from white noise. 
3. To provide for a gust alleviation system which uses readily obtainable 
infomation, and still produce satisfactory performance, the system 


presented in Fig. 4 is proposed. 
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Table TI 
TWO-PARAMETER LONGITUDINAL AUTOPILOT 








Cruising Flight An? 
ET K, l, n^ ДА pO 
Rad/Ft/Airs?/g Rad/Rad/Airs Ft (Ft/Sec*/g)? (Rad)? 
. 36 675 500 . 0466 0022285 0027 
29b 475 1000 0277 ΠΡ, ΓΤ 
5516 400 2000 BRESCH «000818 .41 
36 375 3000 ον ο «000501 «24 
сев 350 4000 ‚00790 .000935  .HH 
. 250 5000 .00637 .000348 Lu 
196 350 6000 . 00531 «000290 .44 
Landing Approach 
An? 
K K^ L, n2 n2 ne SUR 
. 74 2 500 . 0384 . 0138 ΟΙ 
. 711 200 1000 .0209 0069077 230 
WA 175 2000 01.09 „00351 ΠΠ 
. 71 ΠΠ 9000 .00736 00234 87 
«7η 1875 4000 «00556 .00176 «398 
.7Ч 175 5000 ‚00447 ‚00141  .38 
.74 17/5 6000 „00373 OO 239 
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Table IV 
THREE PARAMETER LONGITUDINAL CONTROLLER 
Cruising Flight 


Ki 


K2 


o = 10 Ft/Sec 


K3 


24 


2 


2 


Rad/Rad  Rad/Rad  Rad/Rad (Ft /Sec?/g) ? (вад) 2 
Alrs 
1.59 688.0 Zen .0461 .0001 
1.60 688.0 -2.57 .0324 .0001 
02795 688.0 zu OE .0207 * 
0.568 688.0 το ος .0150 - 
0.452 688.0 EE . 0120 * 
0.393 688.0 2096 «0100 * 
0.447 688.0 — .0085 = 
о = 10 Ft/Sec 
* — Value less than .0001 
Landing Approach 
Kl K2 K3 n? т? 
Rad/Rad  Rad/Rad  Rad/Rad  (Ft/Sec*/g)* (Rad) 2 
Airs 
0.651 400.0 —1.00 .0356 .0011 
0.785 400.0 -1.13 20219 «0008 
255 400.0 U96 ‚0124 «0004 
ο τος 400.0 ο 930 ‚0086 ‚0003 
0.165 400.0 η «0066 «0003 
0.098 400.0 -0.974 .0054 .0002 
0.085 400.0 NOS «0045 «0002 
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FIGURE 1 


Iwo Parameter Gust Alleviation System 





FIGURE 2 


Optimum Gust Alleviation System 
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FICURE 3 


Three Parameter Gust Alleviation System 
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FIGURE 4 
Proposed Gust Alleviation System 
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FIGURE 7 
Three Parameter Gust Alleviation System 
ROOT LOCUS 
CRUISING FLIGHT 
Kl = το... 


Li = 1000 ft. 
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FIGURE 6 
Three Parameter Gust Alleviation System 
ROOT LOCUS 
LANDING APPROACH 


AE το Кз = 1,15 


Lo 1000 ft. 
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APPENDIX A 


Aircraft Characteristics 


The aircraft utilized in this report was the large, rigid, jet 


transport used by Etkin to illustrate application of theory he presented 


in Chapter 6 of Ref. 4. 


The cruising flight data were obtained from 


Ref. 4 and the landing approach data were obtained from Ref. 7. 


Parameter Cruising Flight Landing Approach 
W 1000,000 lbs. 1000,000 lbs. 
S πο sa τε. 1667 sq. fE- 
A 7 7 
Zt. 307000 TE. ο ες 
U. 500 mph. or 733 fps. 200 mph. or 294 fps. 
р 0.000889 slugs/ft.? 0.002377 slugs/ft.? 
c/2 7 7 Ir: Wed oes 
u = m/psc/2 272 102 
= c/2U, 0.0105 sed. 0.0262 sec. 
La B/psc/2 1900 710 
C. 0425 0.59 

O 

0.0188 0.0295 

O 
с, -4.9 -4.8 

a 
C, 0 0 

à 

a 
C -4.2 -4.2 

a 
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Aircraft Characteristics (Continued) 


Parameter Cruising Flight Landing Approach 
С -22.9 -22 ‚9 
m 

d 
8 0 0 

2 

d 
С, -0.24 -0.24 

n 
€ Um 二 ον 
Bn 
C. 0 0 

m 

n 
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APPENDIX B 
Random Gust Fields and Their Use in Autopilot Design 


Aircraft controllers designed to behave well when disturbed by 
sinusoidal, step or impulsive gust disturbances often perform poorly 
when subjected to actual atmospheric conditions. This has led to 
seeking better methods of defining the turbulence or gust fields. It 
has been discovered that if the atmospheric disturbance was treated as 
a random process and the ejrcraft controller designed using statistical 
concepts many of the above difficulties could be eliminated. Controllers 
could then be analytically designed, and built, which performed well 
when subjected to actual atmospheric conditions. 

Statistical design techniques brought with them many new terms and 
concepts not used in the era of the so-called classical autopilot design. 
Terms such as stationary, homogeneous, isotropic, en turbulence, 
mean-square error, correlation functions and power spectral densities 
are now included in the analytical designer's vocabulary. This writer 
will not attempt to provide any mathematical derivations of these 
concepts and will not even discuss them at length. The interested 
Yeader is referred to the bibliography found in Ref. 4, pp. 339-340. 

A very readable review of this subject is found in Ref. 1, Chapter 9. 
In order to utilize statistical techniques in gust alleviation 
design an expression for the power spectral density of the atmospheric 
turbulence is needed. Such expressions have been developed which show 

good agreement with spectra measured in cumulus clouds, thunderstorms 


and clear air. The spectra used for the random vertical gust field of 
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this report is found im Ref. 3, given by: 


il: 
= 2 
acy (2) те (L 9? (1) 


where BS is the mean square value of the vertical gust velocity. L 15 
the characteristic gust length, and © is the spatial frequency. 


Equation (1) can be written in non-dimensional form as: 


2 
O 
wa? 1 (2) 


ου. ος (274125 


z 


Although this expression is not as comprehensive as many expressions 
found for Sg it is used because it is easy to manipulate and it agrees 
fairly well with measured atmospheric power spectra. A form of 
Equation (1) was originally suggested in Ref. 8. d has been found to 
vary from 0 to 40 feet per second and Lo from 500 to 6000 feet. 

The power spectral density of any of the aircraft longitudinal 


motion variables such as angle of attack (a) can be obtained: 
ο D 
$, Loi ο | ва | (3) 


Неге = (ju) represents an aerodynamic transfer function, with s = Ju, 


g 
and relates angle of attack and gust vertical velocity. 
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In addition, mean square values of longitudinal output variables 


can be obtained as: 


1 f oe 
П 8 
EMO (4) 
m Ec ο 
στ j tg ا‎ Aw (5) 


Equations (3) and (5) provide the necessary tools for simple analysis 


and parameter optimization studies. 
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APPENDIX C 
Determination of Transfer Functions 
2 (ju) and 2 (ju) 
"o "g 
I. Longitudinal Equations of Motion. 
The aircraft longitudinal equations below were obtained from Ref. 4, 
page 133: 


X force: = = 72 = 
(2uD as Cy 2а Gef 十 CT 6 0 


“Forces (20, - C, ju + (2uD - Cz D 一 Cz Ja = [ (2u + C, )D 


O u 
= Сна tane 18 ~ Сап = 0 S 


Moment: u ο eu rad opc u) 


= Der )n = 
ng, P+ Gn, 
q - ре = 0 t* = AUN б = с/2 


Substituting a at os q = q + q>) and p=p=p > 
Neglecting changes in forward velocity and considering the X-force 
equation to be identically satisfied yields Equations (1) and (2) below. 
The equations are shown as they appear after the Laplace transform has 
been taken with zero initial condition assumed: 


(24S - C, s ~ C, J)a ~ (2u + C, )q - ee IC, s ΙΕ slag (1) 
α α d α α q ^ 


Е GEN )a + (ips eN la - Ic 5 9 In (2) 
0 0 d n n 


= [С S*C -C slo 
We n, mq 可 
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mete: 4а) ата = рө 


k Alu )dt Um 
πα. 
dy at dx de 


^ 


dg 7 ي‎ 


(р) The = terms above are included on the right side of Equations 
(1) and (2) to indicate their use as forcing functions. 


II. Acceleration Equation. 
The equation for the linear acceleration of the aircraft center of 
mass in the z direction is obtained from: 


а = W + PV - QU 
2 
by neglecting velocities in the y direction, assuming small disturbances, 


and using u = u sina Equation (3) is obtained: 


UÚ ^ 
ac. = cx (Da - q) (3) 


Now one can define the normal acceleration factor: 
-A 


^ c, 2α 5 х 
n (t) = d тт pes (q = Da) 
or in Laplace form 
η. 
S) = WS (q 7 Sa) 4 
DEN = (4) 
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III. Mathematical Model of Longitudinal Control System. 
The equation for the longitudinal control system was obtained from 


the block diagram of the system and is given by: 
Kla + K2q + (K3 - TS - 1)n = 0 (5) 


Note that the elevator-serro combination is approximated by a first 
order system with a time constant (T) of 1/10 second. T in the 


non-dimensionalized form is given by: 


T* = T/t* = 





IV. Simultaneous Solution. 
From Equations (1), (2), (4), and (5) a system of equations is 


Obtained of the fon: 


Aji + Aj od + Άι ηη + A n = e | 
O ο (6) 
A319 + А. + А-1 + Asan = a 
Алу + A23 + A43۳ + А4 П = Όλος 
where 
= S 一 S - = 

Bo: c ETE μα 
A3 = ¬ Cz, Άρη 7 "n; 7 On, 
ποσα E E 
b, = C SFC =C S Ὁ πο S+C -C S 

Z. Z 2 m, m 

α α q à α q 
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The transfer function p(s) de, (S) IE 


Eng 


al 
I 


A 


n 
öl 41 


s=jw A] 


ud 


En 


CH 


an 2 jm 

КО 
K2 А, 20 ас 
K3 - T*s - 1 Aya = 0 
0 Ayy = “1 

by = 0 
is given by: 
S=Ju 

Aa Aus Du 
һә Aus b 
Ау Аз By (7) 
һм» А, bA 
Ἢ... 
Ban Aan Aa 
EEN 
SDN EU 


nie. (5) | is obtained in a similar manner. 


5=)ш 
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APPENDIX D 


Computer Program 


I. The computer program was written to provide a solution to a performance 


index given by: 


со со 


-x ل‎ з AA Dk Go [38 (η) 


со со 


In 


where 


I i E c?c l 
τν е a 
wo 


EIS 


2 (7ш) ` WAN jo) ES (3) 
0 overs 3 + x5 (jw)? + x6(jw) + x 


g 


zl(ju)? 4 x2(ju)? + z3(3u) - z4 и) 
x4 (30)? + x5(30)? + x6(30) + x7 i 





[Gu |? = 
g 


In Equations (3) and (4), xl through x7 and zl through 24 contain 
the variables Kl, K2, and K3 and coefficients of these variables. The 
coefficients are functions of the aircraft stability derivatives and 
the gust. characteristic length. 

It can be seen that the performance index can be written in the form: 

LE Cljw)C(-jo) 


Iz AÑ da 


n 2T" J d(jo)d(-jo) 


(5) 
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This integral is evaluated for n = 1 through n = 10 in Ref. 5. For the 


performance index of this problem 1, = ر1‎ where 


2 2 2 DM 2 2 
3 (And al dp) + (c5 2c)c.)d d.d +(c] 2c ἃ Heg e (daa qui qo) 


Ἢ Baron 02194 00304190 01941900380) e 
a 9 2 
248, (-d.d5s—dyd +A dd) 


c(jw) + c3 (jw) > + co (ju)? + cy (jw) + cp 
d(jw) = dy (jw) * + d3 (jw) SS d» (jo) ? ΟΠΠ da 


The I, equation was programmed and solved on the digital camputer. 


4 
The program was written so that the type aircraft could be changed, i.e. 
all stability derivatives were retained in the computation of xl through 
x7 and zl through z4. ‘The gust characteristic length (CL) was set at 

the beginning of the program. To facilitate handling of the variables 

Kl, K2, and K3 their coefficients were stored in 3 dimensional arrays 
according to the exponents of the variables. For example, the coefficient 
of K12K2K33 would be stored in an array at location E, 2,4). The value 
stored in this location was printed out as PI(2,1,3) = (value stored). 

The arrays were manipulated to solve for the performance index in terms 

of Kl, K2, and K3. 

A minimum value of PI(Kl, K2, K3) was determined by a steepest 
descent method. A starting point for the minimization procedure was 
obtained by computing values of the performance index as Kl, K2, and K3 
varied over a large range. The point where the performance index was 
a minimum was then used as a starting point for the steepest descent 
method. Several other starting points were utilized to help ensure that 


a proper minimum was reached. 
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IT. Flow Chart 








SEP CL, Kl, K2, AND K3 
ARG(1) = Kl; ARG(2) = K2; ARG(3) 






COMPUTE AND SET VALUES OF d AND 
C TO SOLVE THE INTEGRAL 


SEE PAGE 41 OF THIS APPENDIX 
MEANING OF d AND d 










CALL SUBROUTINES TO CALCULATE 
PERFORMANCE INDEX, SEE SECTION III 
OF THIS APPENDIX FOR SUBROUTINE 
EXPLANATION 


WRITE OUT PERFORMANCE INDEX NUMERATOR 
AND DENOMINATOR IN THE FORM: 
Pu 
WHERE I,J, AND K ARE THE EXPONENTS 
OF Kl, K2, AND K3 RESPECTIVALLY. 
PI(2, 1, 3) = .06 WOULD BE .06K1%K2K3? 



















COMPUTE THE PARTIALS OF THE PERFORMANCE 
INDEX NUMERATOR AND DENOMINATOR 
WITH RESPECT TO Kl, K2, AND K3 


CALL SUBROUTINE FUNCT WHICH DETERMINES 
THE VALUE OF THE PERFORMANCE INDEX 
AND COMPUTES THE GRADIENTS IN THE 
K1, K2, AND K3 DIRECTIONS FOR A 
GIVEN Kl, K2, AND K3. THIS 
SUBROUTINE USES THE PARTIALS 
COMPUTED IMMEDIATELY ABOVE. 
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PRINI PERFORMANCE INDEX, 
Kl, K2, K3, AND GRADIENTS 
IN Kl, K2, AND K3 DIRECTIONS. 


s < YES 
ALL GRADIENTS LESS THAN 107% > 
STEP SIZE LESS THAN 10-!6 FT 
NO 













STEP IN DIRECTION O 
STEEPEST DESCENT 


YES = 
ANYMORE CONSTANTS 
GL, Kl, K2, AND 7А 


NO 


» 
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IT. 


Subroutines 


ZERO3 
ZERO4 


MULT1 


MULT2 


MULT3 


SUBTI 


SUBT2 


STUFEP 


- Places zeros in a 3x3x3 matrix. 
- Places zeros in a 4x4x4 matrix. 


- Multiplies a first order polynomial by a first order 


polynomial and stores answer in a 3x3x3 matrix. 


- Determines the product of a first order polynamial 


times a first order polynamial times 2.0 and stores 
answer in a 3x3x3 matrix. 


- Multiplies a first order polynomial by a second order 


polynomial and stores answer in a 4x4x4 matrix. 


Subtracts a third order polynomial from a third order 
polynamial and stores answer in a 4x4x4 matrix. 


Subtracts a second order polynomial from a second 
order polynomial and stores answer in a 4x4x4 matrix. 


Determines the product of a second order polynomial 
times a third order polyncmial times a constant and 
stores the answer in a 6x6x6 matrix. Separate 6x6x6 
matrices hold the numerator and the denominator of 
the performance index. 


Takes the performance index numerator and denominator 
arrays, the specific values of Kl, K2, and K3, and the 
arrays holding the partials of the performance index 
numerator and denominator and calculates the value of 
the performance index and the gradients in the Kl, K2 
and K3 directions. 
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